Matrix Product States and the Fractional Quantum Hall Effect 
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We obtain an exact matrix-product-state (MPS) representation of a large series of fractional 
quantum Hall (FQH) states in various geometries of genus 0. The states in question include all paired 
k = 2 Jack polynomials, such as the Moore-Read and Gaffnian states, as well as the Read-Rezayi 
k = 3 state. We also outline the procedures through which the MPS of other model FQH states can 
be obtained, provided their wavefunction can be written as a correlator in a 1 + 1 conformal field 
theory (CFT). The auxiliary Hilbert space of the MPS, which gives the counting of the entanglement 
spectrum, is then simply the Hilbert space of the underlying CFT. This formalism enlightens the 
link between entanglement spectrum and edge modes. Properties of model wavefunctions such as 
the thin-torus root partitions and squeezing are recast in the MPS form, and numerical benchmarks 
for the accuracy of the new MPS prescription in various geometries are provided. 



The understanding and simulation of quantum 
many-body states in one space dimension has experi- 
enced revolutionary progress with the advent of the 
density matrix renormalization group pQ. In modern 
language, this method can be viewed as a variational 
optimization over the set of matrix product states 
(MPS) |3]. Indeed, gapped one-dimensional sys- 
tems (which generally have low entanglement) can be 
very efficiently simulated by expressing the weights 
of many-body non-interacting states in an interacting 
wavefunction as products of finite-dimensional matri- 
ces B[i] mi associated with each occupied (m, = 1) 
or unoccupied (m, = 0) site i (or spin): = 
E {ms} P L B[l] m iB[2} m * . . . B[N s ] m «sP R |m 1 ,m 2 , . . . m Ns ) 
Pl,Pr are projectors into the state at the end of the 
chain (absent for periodic boundary conditions or the 
infinite chain). As long as the "bond dimension" \ 01 
the matrix B[i] is less than 2 l , this provides a more 
economical representation of the state. Generic 1-D 
gapped systems can be approximated by finite \ [I]- 

Due to their perimeter law entanglement, 2-D systems 
(such as the fractional quantum Hall effect) are harder 
to simulate by MPS [5--8J. In a recent paper [9], ex- 
ploiting the fact that continuum model FQH states can 
be written as correlators of primary fields in conformal 
field theories (CFTs), an MPS expression was obtained 
for continuum Laughlin [TO] and Moore-Read [TTj states 
on infinite cylinders. The bond dimension \ grows with 
the number of particles but scales with the circumfer- 
ence L of the cylinder rather than its area. Approximate 
expressions can be obtained by truncating \ of the ex- 
act MPS. A key ingredient of Ref. [5] is the expansion of 
operators in a free basis (boson for Laughlin and Majo- 
rana plus boson for Moore- Read) , which cannot be easily 
implemented in the more complicated, interacting bases 
of other FQH states such as the Read-Rezayi series [12]. 
The exciting possibility is that, if all model FQH states 
could be written in MPS form, current numerical bar- 



riers could be broken and properties such as correlation 
functions would be computable for large sizes. This is 
supported by the continuous MPS proposed in Ref. |13) . 

In this paper we provide a generic prescription that en- 
ables us to obtain the unnormalized (thin annulus) MPS 
form of a model FQH state which is the correlator of a 
primary field in a CFT. We explicitly construct MPS for 
the (k,r) — (2,r) paired Jack polynomial states [Til [TO] 
(r = 2,3 being the Moore- Read and Gaffnian [TO] wave- 
functions), as well as for (fc,r) = (3,2) corresponding 
to the Read-Rezayi Z% wavefunction [T5]. Several key 
ingredients and subtleties, such as the presence of non- 
orthonormal bases, null vectors, and intricate operator 
commutation relations, are discussed. We then show how 
to extend the MPS description to different manifolds such 
as the cylinder, sphere, and plane, and how several known 
properties of the CFT wavefunctions such as squeezing 
arise naturally in this description. We then generate sev- 
eral (the (k,r) = (2,2), (2,3), (2,6), (3,2)) of these states 
numerically to verify our MPS, and provide numerical 
benchmarks to attest the accuracy of the MPS on the 
cylinder and sphere. 

A large class of FQH ground-states are described by 
the many-point correlation function of an electron oper- 
ator field V{z) in a chiral 1 + 1 CFT [TT]: 

c\m x (z l , . . . ,z N J (1) 

A 

where m\ are monomials (or Slaters for fermions) of an- 
gular momentum A = (Ai, . . . , Ajv c ). N e is the number of 
electrons, and the filling fraction is v — 1/q (note that q 
will not always be integer). The state (N ey /q\ describes 
the background charge at infinity. The coefficient c\ can 
be obtained by contour integrals. Upon inserting a com- 
plete basis of states in the l.h.s. of we get 

^En^ ( 2 ) 

{aj}j=l J Z 3 



2 



0, N ey /q, re- 
site (Landau 
MPS I*) = 



mi-.- m N<t 



with 



The J7(l) charge of \ao) , |o;Ar c ) is 
spectively. This is an infinite, 
level orbital momentum) dependent 

the matrices for an orbital j being, in the limit of an 
annulus with a very large radius (the so-called "con- 
formal limit"), (a'\B°\j]\a) = 6 a ,, a and (a'\B 1 [j]\a) = 
<5a ,.A a +h+j(ct'\V(l)\a}. h is the conformal dimension 
of V(z), and higher occupation number (of occupation 

m) matrices are simply B m [j] — [B 1 ^} ) /ym!. 

To obtain a site-independent MPS, we need to spread 
the background charge uniformly over the droplet. We 
make explicit the dependence on the U(l) charge by writ- 
ing states | a) = |Q)<8>|a), where Q is the U(l) charge and 
a encodes the rest (descendant, neutral sector). As the 
matrix element (a' 1^(1)] a) does not depend on charge Q, 
we are free to modify the distribution of the background 
charge. Spreading uniformly the background charge 
amounts to an insertion of a U(l) background charge 
— 1/y/q between each orbital. This yields a site inde- 
pendent MPS with B m = e -</2v^oym e -i/2Vg<P0 / v / TO T j 
where (pa is the U(l) boson zero mode. 
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where a and a' are the basis of descendants in the CFT 
(not free fermions as in Ref. Our expression further 
differs from the one in Ref. |9] by time-evolution terms 
U(St) = exp(— StLq), which give the cylinder normaliza- 
tion. Conformal invariance yields 



(a'\V(z)\a) =z L 



- h (a'\V(l)\a), (4) 



where h, A a are the conformal dimensions of the primary 
field V(z) and of the descendent \a). The matrix ele- 
ments of Vq are then related to the CFT 3-point func- 
tion (a'\V \a) = 6 Aa ,Ac,+h(a'\V{l)\a}. The "electron 
operator" is a primary field of the tensor product form 
V(z) = $(z)(8s : e^vW :, where $(z) lives in the so- 
called neutral \a) conformal field theory CFT„ factorized 
from the U(l) sector \b). In this basis \a) — \a) ® \b) 
the 3-point function factorizes as (a';b'\V(l)\a;b) = 
5 Aa ,+A b ,,Aa+A b +h(a?\Hl)\a) (V\ : e^W : |6> where 
h = q 2 /2 + /i$ is the conformal dimension of V(z). Note 
the delta function in the total conformal dimension of the 
field and not separately in the neutral and U(l) parts. In 
the following we explain how to obtain the neutral, inter- 
acting CFT matrix elements for the case of (k, r) Jacks 
for k = 2 and k = 3. 

First, we re-evaluate the U(l) matrix elements in a way 
easily generalizable to non-free field CFT and in a basis 
where they are real: 



ip(w) =<po~ ia o log(w) + i y~] -a„w n , 



(5) 



where a n are the bosonic modes obeying the Heisenberg 
algebra [a n ,a m ] = n8 n +m,o- is the zero mode of the 
conserved current and measures the U(l) charge, while 
ipo is its canonical conjugate ([(po, do] = i). Primary fields 
are the vertex operators Vp(z) —: e l ^ v ^ : with conformal 
dimension /3 2 /2. The corresponding highest weight state 
1/3) = V^(0)|0), which is annihilated by all a„>o, has 
U(l) charge f3. Since we defined the U(l) charge to be 
a , the electric charge is ^Qo- Descendants are obtained 

by acting on |/3) with the lowering operators a\ = a_„, 
n > 0. They are labeled by a partition /i = {/ij}, with 

Ml > A*2 > ' ' ' > > 

\Q,H) = Y[a^\Q), a \Q,ri=Q\Q,») (6) 



For multiplicities of element j, m,j — mj(fj,) the norm is: 
(Q,M|<2V) = ^<W<W> z li = Y[j m *m j \ (7) 

3 

The matrix elements of the primary field between the 
normalized basis of descendants can be computed easily 



using the recurrence [a m , : e 
They are of the form 



a if3ip(z) . p z m 



(8) 

which is consistent with Eq. Q and with charge conserva- 
tion; f3Q comes from the difference in conformal dimen- 
sions (Q' 2 — Q 2 )/2. With rrij = rrij(fi) and m'j = rnj(^i') 
being the multiplicities of j in fi and fi' , one finds: 



A„ 



j>l r =0 s=0 



m'j \rrij ! 



(rrij - s)\ 
(9) 



These matrix elements are real (useful for numerics) with 
the symmetry A^ifi) = A^ l(tl (- ( 8). 

We now move to the neutral part, starting with the 
case of (k,r) = (2,r) Jacks. The CFT can be factor- 
ized as a U(l) free boson times a neutral minimal model 
M(3, 2 + r) |17j . The underlying symmetry of this mini- 
mal model is the Virasoro algebra 

[L n , L m ] ={n- m)L n+m + ^-n(n 2 - l)<5 n+m! o (10) 
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with central charge c = 1 — Q(g — l) 2 /g, g 
The electron operator is V(z) = : e lV2h "- +m ^ z l 

^(z) = $(!| 2 )(2) is a primary field in the neutral CFT, 
with conformal dimension h^,, kept generic. Like in 
the U(l) case, the Hilbert space of the neutral CFT is 
made of the primary fields |A), eigenstates (with eigen- 
value A) of Lq and annihilated by lowering operators 
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L n , n > 0, and their descendants, indexed by a partition 
A = (Ai, A2, • • • , A„): 



|A,A) = i_ Al i-A 2 ---i-AjA), 



(11) 



also eigenstates of L with eigenvalues A + |A|, where 
|A| = Using the Virasoro algebra (10) and 

L\ = L- n , we can compute any overlap between descen- 
dants. While descendants with different |A| are clearly 
orthogonal (having different Lq eigenvalues), a major dif- 



ference from the 17(1) case is that the descendants (11) 
are generically independent but not orthogonal. Hence, 
at each level |A|, we have to numerically build an or- 
thonormal basis. 

Another issue is that for non unitary CFTs (such as 
the one underlying the Gaffnian) some descendants have 
a negative norm. These states have to be included, and 
the sign of their norm can be handled by an extra diag- 
onal matrix D acting on the right of the MPS matrices : 
B m ^ B m D. 

A final major difference from the C/(l) case is the pres- 
ence of null vectors - states of vanishing norm under the 
scalar product defined by L\ — L- n . This is a reflection 
of the fact that for special values of A (which include 



all the interesting cases), some states in (11) are not in- 
dependent. CFT characters count the number of inde- 
pendent descendants at each level, from which one can 
deduce the number of null- vectors. Our numerical pro- 
cedure for computing overlaps reproduces this counting. 
At each level we need to detect and drop all null-vectors 
before performing the Gram-Schmidt process. 

To compute matrix elements (A', A'|^I / (1)|A, A) be- 
tween descendants, it is convenient to work in the over- 
complete (due to null vectors) "basis" (11) and then 
transform back to the orthonormal basis. The level-0 
matrix element (A'|^(1)|A) is simply the OPE structure 
constant D^>^^,A, known in the closed form for minimal 
models and, gives an overall pre-factor which can be ig- 
nored. Others can be computed using a similar method 
as for the U(l) CFT; for all m £ Z, 



[L m -L ,*W(l)]=mh& h \l), 



(12) 



where is any primary field with a generic conformal 
dimension h. Any matrix element can in principle be 
computed exactly using this method but to the best of 
our knowledge there is no analytical closed formula. 

We are now in the position to write down the MPS ma- 
trices for (A; = 2, r) Jack states. The state | A, A; Q, fi) = 
I A, A) <E> \Q, fj) of an over-complete "basis" of states in the 
CFT for k — 2 Jack states has level ("momentum") P = 
|A| + |/x|, which serves as the truncation parameter for the 
MPS. The matrix elements (A', A'; Q', ^'|B m | A, A; Q, fi) 
for m = 0, 1 are given by: 



B 1 '■ ^A'+|A'|+Im'I+V90.A+|A|+|/*|+/»*+i/ 2 
x{A',X'Ml)\A,X)S Q , <Q+ ^_ 1/ ^A fi ,^. 



(14) 



B a 



W'*Q',Q-^A,A'<A',A'|A,A) 



(13) 



(A', A'| 1 I'(1)|A, A) can be computed using the neutral 
CFT, then changed to an orthonormal basis. A^i^ is 
given m ^ for $ = y/q. The values of A, A' are fixed 
by the fusion rules of the electron operator in the neutral 
sector. For the (k, r) — (2, r) Jacks, these are $X$ = 1, 
f x 1 = which gives rise to two neutral sectors \x) 
x = 0, 1 with conformal dimension A = xh^, and 

(x',\'\*(l)\x,\) = 6 X+X , tl (x',\'\*(l)\x,\). (15) 

The charge U(l) sector is specified by the possible values 
of Q, Q'. As q = 2h\s + m = r/2 + m, m € Z, Q, Q' are 
always shifted by integer (for r even) or half integer (for r 
odd) multiples of 1/y/q, so Q = + Q , n € Z, |Z for r 
even, odd. The value Qo £ K can be determined from the 
requirement in (III]) that A' + \\'\ + \n'\ + y /qQ = A + |A| + 
|/i| + /i3< + l/2, where A = xh^, A' = x'h<$, and x' = 1—x. 
This yields the constraint n + ^/qQo = (xr + l)/2 mod 1. 
For r odd, this also implies an extra constraint between 
the neutral sector x and the charge sector 

1 _ 1 — x 

V^QeZ + -, reven; v/gQeZ + — , r odd. (16) 

We have implemented the above MPS numerically 
and verified that it exactly reproduces all (fc,r) = 
(2, 2), (2, 3), (2, 6) Jack states. As a short exercise, it is in- 
structive to show that P = MPS recovers the thin torus 
limit [TS] (root partition) of these Jacks. Denoting N — 
y/qQ, we find B°, N ,. x N oc 6 x > tX 5jv',jv-i, b I>,n'-,x,n k 
&x'+x,i ^x,(2N-i)/r^N',N+g-i from which we see that B 1 
can act only on two states : \x = 0, N = 1/2) and 
\x = 1,N = (r + l)/2). Starting with the state \x = 
0, N = 1/2) we get the correct root partition for the 
(fc, r) = (2, r) Jack polynomial multiplied by m Jastrow 
factors: (. . . r+m - 1 10 m - 1 10 r+m - 1 10 m - 1 l). The topo- 
logical sector (responsible for the ground-state degener- 
acy) can be fixed by choosing matrix elements of the 
product B mi . . . B mN *> between different primary fields. 
One can describe quasi-hole states by inserting quasi-hole 
matrices in the MPS. Or alternatively, edge states are ob- 
tained by choosing matrix elements involving descendant 
states instead of primaires. This means that the MPS for- 
malism establishes a mapping between edge states and 
the auxiliary space, which in turns controls the entangle- 
ment spectrum. In particular the MPS makes transpar- 
ent the counting of the orbital entanglement spectrum 
for such states [!§] . 

We now move to the Read-Rezayi Z 3 state, exempli- 
fied by the k = 3 Jack polynomial. For k > 2 Jack 
states, the only known approach is to deal with a CFT 
with an enlarged algebra (the so-called Wfc algebra) [see 
Eqs. (44), (45), and (46) of Ref. 20^ which includes a W 
current of spin 3. This generic approach, which applies 
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to all (k, r) = (3,r) Jack states becomes inefficient for 
the Z 3 RR state due to the appearance of an extremely 
large number of null vectors (at each level of truncation). 
The underlying CFT for the k — 3 RR state is known 
to be equivalent to the minimal model M(5, 6), the field 
^i(z) becoming the primary field ^^^(z). This alterna- 
tive approach provides a basis in which matrix elements 
involve only Virasoro modes L n . The electron opera- 
tor is V{z) =: eVwM . ^^(z) with q = 2/3 + m. 
The neutral CFT field $(3|i) can be split into two chiral 
fields ^i(z), 4'_x(z) with conformal dimension = 2/3. 
Their fusion rules in the W3 framework are f 1 x = 
1 J r _i,l'i x = 1. While in the W3 algebra language 
\W) = ^/3/cW-3 |0) is a descendant of the identity, it 
is primary (with conformal dimension 3) with respect to 
the Virasoro algebra: L n \ W) = 0, n > 0. Accordingly, 
in the minimal model M (5, 6) framework one has to work 
with the fusion rule x W_i = 1 + W. 

The Z3 parafermions have three sectors correspond- 
ing to the Z 3 charge of the field x = 0,±1. Work- 
ing in the Virasoro algebra, the x = sector contains 
two primaries |0) , \ W) as well as their descendants ob- 
tained just like above by the action of the Virasoro gen- 
erators P_„, whereas x = ±1 are made of ^±1) an d 
their descendants. The matrix element between descen- 
dants (A', A'|\E f i(l)| A, A) vanishes unless x' = x + 1 mod 
3. The matrix elements we need are (\J>i, A'|\I/i(l)|0, A), 
($1, A'|*i(l)|W,A), (*_i,A'|*i(l)|*i, A), all other be- 
ing obtained from the above by charge conjugation 
= (C(a)\^i(l)\C(a')), where charge con- 
jugation interchanges ^±1), leaves |0) invariant and 
flips the sign of \W). 



we can coin- 



Using Eq.|12j), 
pute these matrix elements up to one coefficient, namely 
*i(l) \W) I *i(l) |0). But this is simply an 
OPE structure constant, which is found to be \/26/9. 

Once the matrix elements (real, with this normal- 
ization) between descendants are known, it is easy 
to find the explicit form of the MPS B matri- 
ces for the RR state, (A', A'; Q' , fJ,'\B°\ A, A; Q, /jt) and 
(X,\';Q',H , \B 1 \A,\iQ,n): 



B° 



'W'^Q'.Q- 



B 1 



^A' + |A'| + | M '| + x /gQ.A+|A| + | AI | + h 4 , + l/2 



(A',A'|*i(l)|A,A) Sq, 



with = 2/3. is given in ^ for /3 



1/V? <A',A'|A,A) (17) 

+\\\+\n\+hv+l/2 

(18) 

Vq- 

A = {Xi} is a partition of descendants of all the four 
primary fields (W included) in the theory. As before, 
P = |A| + \fi\ is the truncation parameter for the MPS. 
As an exercise, we obtain the P — thin torus limit of 
the MPS. The W field will play no role for P < 2, and 
there is a unique primary field per sector (for P > 2 
there are two primary fields in the sector x = : |0) and 
\W)). At level P = we find B° = S x , , X S N , , N ^, B 1 = 

Sx',x+l$N>,N+q-lS2\x'\/3+N,2(\x\ + l)/3+l/2 w hh 1 = 2/3 + 




Laughlinv=l/3 — •— 
Moore-Read (k=2 r r=iz) — •— 

Gaffnian (k=2,r=3) «--- 

Read-Rezayi k-3 — *— 
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FIG. 1. (Color online). Growth of the auxiliary Hilbert space 
and convergence criteria for different types of boundary con- 
ditions. Upper panel: Dimension of the auxiliary space for the 
Laughlin, Moore- Read, Gaffnian and Read-Rezayi Zz state as 
a function of P ma x- Lower panel : Quantum distance between 
the exact Read-Rezayi Z3 and the truncated MPS state on 
the sphere (b) and on the cylinder with aspect ratio one (c) as 
a function of Pmax- Note that for N — 24 and P max = 10, the 
MPS only gives 40% of the 6.2 x 10 7 components. Smaller 
quantum distances are obtained for Laughlin, Moore-Read 
and Gaffnian (not shown). 



m. Acting on the state \x = 0, N = 1/2) we recover 
the root partition . . . 10 m - 1 10 m - 1 10 m+1 10 m - 1 10 m - 1 l, 
which is the thin torus limit of the RR state multiplied 
by m Jastrow factors. 

The obtained MPS description of the Jacks (un- 
normalized wavefunctions on the annulus) is trivially 
transmuted to other geometries: for the sphere and the 
infinite plane, we obtain a site dependent MPS with 
.B mi [i] = B mi /N(i), where iVj is the norm of the sin- 
gle particle orbital z l in the respective geometries. For 
the cylinder, a site independent MPS is possible by intro- 
ducing a time-evolution e ~ 27rL o/L^ w j^ n £ denoting the 
circumference of the cylinder. The well-known squeez- 
ing properties of FQH states easily follow from the MPS 
description. 

So far the MPS description we have provided is exact. 
For numerical purposes of approximating a given state, 
we introduce the truncation level P ma x, which is the max- 
imum allowed value of P = |A| + P m ax = gives 
the root partition (thin-torus limit) of the Jack polyno- 
mials and amounts to dropping all descendants in the 
CFT Hilbert space. P m ax — 1 gives the correct weights 
for all configurations obtained through a single squeez- 
ing of the root partition. Generically, the truncation at 
P max amounts to restricting the number of any squeezing 
steps from the root partition to P max - This is also the 
momentum quantum number labeling the entanglement 
spectrum levels [5T] . Due to the shape of the orbital spec- 
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trum, the truncation to a certain P max is expected to be 
equivalent to keeping the states with the highest Schmidt 
weight in the groundstate. This is also related (though 
not equivalent) to expansions around the thin cylinder 
limit [35], where the weights of configurations decrease 
with the amount of squeezings from the root partition 
(equivalent to the exponential decay of correlation func- 
tions in the associated CFT). 

In Fig. [T] we provide numerical benchmarks for the ac- 
curacy of approximating the full Jack states by the MPS 
truncated at level P max . The approximate Laughlin, 
Moore-Read, Gaffnian and Read-Rezayi states have been 
constructed by MPS B matrices whose auxiliary Hilbert 
space dimension grows as shown in Fig. [TJa) . The accu- 
racy of the approximation is quantified by the overlap of 
an MPS state with a full Jack polynomial, and an exam- 
ple for the Read-Rezayi Z3 state is given in Figs. [ljb,c). 
We observe that the approximate MPS state becomes 
an excellent approximation of the exact FQH state for 
relatively low values of P max = 10. Note that the con- 
vergence to the exact state on the sphere [Fig. |T|(b)] is 
strikingly slower than on the cylinder [Fig.JTJc)], making 
this a preferred type of boundary condition for DMRG 
implementations (23J [24] • 

In conclusion, we have provided a method to obtain the 
MPS description of FQH model states given by correla- 
tors of the primary fields in a CFT. We have furthermore 
obtained the exact MPS form of the (fc, r) = (2, r), (3, 2) 
Jack states (including Moore-Read, Gaffnian, and Read- 
Rezayi k = 3 state) , and compared the approximate MPS 
(truncated to a certain P max ) with the exact states. Com- 
paratively small values of P max were found to be suf- 
ficient for obtaining extremely accurate approximations 
of these states on the cylinder, which might have con- 
sequences for the improved DMRG implementations of 
realistic (Coulomb) Hamiltonians. 
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